arXiv:1503.07941vl [physics.atom-ph] 27 Mar 2015 


Frequency Ratio of and Optical Lattice Clocks beyond the SI Limit 


Kazuhiro Yamanaka/’^’^ Noriaki Ohmae/’^’^ Ichiro Ushijima,^’^ Masao Takamoto,^’^ and Hidetoshi Katori^’^’^ 

^Quantum Metrology Laboratory, RIKEN, Wako-shi, Saitama 351-0198, Japan 
‘^Innovative Spaee-Time Project, ERATO, JST, Bunkyo-ku, Tokyo 113-8656, Japan 
^Department of Applied Physics, Graduate School of Engineering, 

The University of Tokyo, Bunkyo-ku, Tokyo 113-8656, Japan 
(Dated: March 30, 2015) 

We report on a frequency ratio measurement of a ^^^Hg-based optical lattice clock referencing a 
^^Sr-based clock. Evaluations of lattice light shift, including atomic-motion-dependent shift, enable 
us to achieve a total systematic uncertainty of 7.2 x 10“^^ for the Hg clock. The frequency ratio is 
measured to be i^Hg/z^Sr = 2.629 314 209 898 909 60(22) with a fractional uncertainty of 8.4 x 10“^^, 
which is smaller than the uncertainty of the realization of the SI second, i.e., the SI limit. 

PACS numbers: 06.30.Ft,37.10.Jk,32.60.+i, 42.62.Eh, 42.62.Fi 


Rapid progress in optical lattice clocks [iHi makes 
them potential candidates for the future redefinition of 
the second and new tools for testing the fundamental 
laws of physics [^, 0- The absolute frequencies of ^^Sr 
and ^^^Yb-based optical lattice clocks have already been 
measured with the uncertainty of the realization of the 
SI second, allowing them to be adopted as the secondary 
representations of the second by Comite International des 
Poids et Mesures (CIPM) Q. Moreover, the systematic 
uncertainties of the ^^Sr clock reach 10“^^ level Hi, 
which is competitive to ion-based clocks j^. Such clocks 
can be used as new references to investigate precise fre¬ 
quency measurements far beyond the SI second. 

Hg is another promising candidate for optical lattice 
clocks [Sliil, because its susceptibility to the blackbody 
radiation (^BR) is an order of magnitude smaller than 
that of Sr [l-[^ and Yb 0- Furthermore, a large nuclear 
charge of Z = 80 makes the Hg clock a sensitive probe for 
testing the constancy of the fine-structure constant 0 . 
So far the absolute frequency z^Hg of the ^ So — ^Po clock 
transition of ^^^Hg has been reported with an uncertainty 
of 5.7 X 10 [ll|, where the SI second is sufficient to 

describe the frequency. 

In this Letter, we report on a frequency measure¬ 
ment of ^^^Hg by referencing to the ^^Sr clock frequency 
^Sr [ll- We have determined the frequency ratio to be 
^Hg/^Sr = 2.629 314 209 898 909 60(22) with a fractional 
uncertainty of 8.4 x 10“^^. This ratio can be converted 
to z/Hg = 1 128 575 290 808 155.4(1.1) Hz via the recom¬ 
mended frequency for the Sr transition given 

by CIPM with uncertainty 1x10“^^. To improve the 
systematic uncertainty of the Hg clock to 7.2 x 10“^^, 
we investigated the lattice light shift taking into account 
the multipolar polarizabilities 0 , which affect the light 
shift more seriously than in Sr and Yb clocks. 

Figure [T](a) shows the experimental setup for the fre¬ 
quency ratio measurements, which consists of a ^^^Hg 
clock, a ^^Sr clock 0, and a frequency link between 
them. The — ^Po transition of ^^^Hg with a nu¬ 
clear spin of / = 1/2 [see Fig. [11(b)] is used as the clock 


transition and is investigated by the following experimen¬ 
tal sequences with a cycle time of 1.5 s. Hg atoms are 
laser-cooled by a vapor-cell type magneto-optical trap¬ 
ping (MOT) on the ^Sq — ^Pi transition at 254 nm with 
a natural linewidth P = 1.3 MHz 0 . Initially, to collect 
atoms, we apply a MOT laser detuning Sumot ~ — 7T 
and an intensity ~ lOmW/cm^ per beam. After 860- 
ms-long atom loading time, we increase the gradient of 
quadrupole magnetic field from I mT/cm to 3 mT/cm 
to compress the atom cloud. We finally tune the laser 
frequency to 6umot ^ —P and reduce the laser inten¬ 
sity to further cool the atoms to maximize the transfer 
efficiency into a lattice trap. 

About 3 % of atoms are loaded into a vertically- 
oriented one-dimensional (ID) optical lattice operated at 
Al ~ 362.6 nm 0 with the maximum trap depth of 
7 /l ~ 65 Fr, where Fr/Zz = h/{2mX‘l) = 7.6 kHz is the 
lattice-photon recoil frequency, h is the Planck constant, 
and m is the mass of ^^^Hg atom. We then temporarily 
decrease the trap depth down to 7 /l ~ 35 Fr to release 
the atoms trapped in high-lying axial vibrational states 
with n > 3. This ID lattice is formed inside a buildup 
cavity with a power enhancement factor of ~ 10, which 
consists of two curved mirrors and a plane folding mirror 
[see Fig. [11(a)]. These cavity mirrors are placed outside 
the vacuum chamber to prevent vacuum-degradation of 
the mirror coating JjJ , while the Brewster windows pro¬ 
vide optical access to the vacuum and selective enhance¬ 
ment of the p-polarized light. Atoms in the lattice are 
spin-polarized by exciting the ^ So — ^Pi transition with 
the circularly polarized light. A bias magnetic field Bbias 
is applied during the spin polarization and clock excita¬ 
tion. 


The ^So — ^Po transition is excited by a clock laser at 
266 nm generated by two-stage frequency doubling of a 
fiber laser at 1063 nm stabilized to an Er-doped fiber op¬ 
tical frequency comb by a linewidth transfer method [l5| . 
The carrier envelope offset frequency /ceo of the comb 
is stabilized by using a self-referencing / — 2/ interfer¬ 
ometer. The repetition rate /rep of the comb is then 




FIG. 1. (a) Experimental setup. Laser-cooled Hg atoms 

are loaded into an optical lattice at the magic wavelength of 
363 nm. The ^Sq — ^Pq clock transition is interrogated by 
a clock laser at 266 nm, which is referenced to the Sr clock 
laser at 698 nm via an optical frequency comb. SHG, second 
harmonic generation; GGD, charge-coupled device; FS, fre¬ 
quency shifter; DDS, direct digital synthesizer, (b) Relevant 
energy levels for ^^^Hg with a nuclear spin of / = 1/2. (c) A 
clock spectrum of the ^ So —^Pq transition of ^®®Hg. 


stabilized by referencing a sub-harmonic (1397 nm) of a 
Sr clock laser at 698 nm, which is prestabilized to a stable 
reference cavity with instability ?^5xl0“^^atr = ls 
and stabilized to the Sr clock transition for r > 10 s 

The clock laser is superimposed on the lattice laser 
with the same polarization. Figure [T]^c) shows a clock 
spectrum with a Fourier-limited Rabi linewidth of yng = 
7.4(4) Hz, corresponding to a Q-factor of z^Hg/7Hg = 
1.5 X 10^^, which is obtained for the clock interrogation 
time Ti = 120 ms. For the data presented below, we op¬ 
erate the Hg clock with ri = 40 — 80 ms so that the fre¬ 
quency stabilization to the atomic transition becomes ro¬ 
bust against the variations of experimental conditions for 
several hours. The atom population Ns in the ^So state is 
determined from the fluorescence by operating the MOT 
for 20 ms. We then optically pump the atoms in the 
^Po state back to the ^ So state by exciting the ^Pq — 
transition at 405 nm to determine Np of atoms in the 
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FIG. 2. Glock frequency shift Az/c(7l, as a function 
of the lattice laser intensity /l (bottom axis). Gorrespond- 
ing lattice trap depth Pl is given in the top axis, (a) Mea¬ 
sured data (crosses) are fitted to Eq. ([T]) shown by curves 
with corresponding colors. Numbers below the curves indi¬ 
cate the lattice laser frequency detuning (5 z/l (in MHz) from 
^ = 826 855 533 MHz. (b) An enlarged view for the vertical 
axis. For (5 z/l = —2 MHz used for the ratio measurements, the 
estimated clock frequency shift and uncertainty are shown by 
a black curve and a gray shaded region, respectively. We op¬ 
erate the lattice at /l ~ 59kW/cm^ or Pl ~ 43Pr as shown 
by a dashed line. 


^Po state. The excited atom fraction is used to 

stabilize the clock laser frequency. We alternately inter¬ 
rogate the two TT transitions mp = ±1/2 ^ mp = ±1/2 
to average out the Ist-order Zeeman shift and the vector 
light shift 0- 

The lattice light shift has been the primary source of 
the uncertainty of ^^^Hg clock [ll|. To evaluate the lat¬ 
tice light shift Az^c(7l, z^l), we measure the intensity- 
dependent clock shift Suc{Il, = Az/c(/l,^l) — 

Az/c(/o 5 ^l) in successive measurement cycles by vary¬ 
ing the lattice laser intensity Ip, while we keep Jq = 
89kW/cm^ constant. Here the intensities Ip are given 
in terms of the peak intensity of a single traveling-wave 
laser that creates the lattice potential depth of Pp. Fig¬ 
ure [2] (a) shows the data taken at 7 different lattice laser 
frequencies z/p, which are stabilized to the optical fre¬ 
quency comb within 0.5 MHz. Each data point is mea¬ 
sured with an uncertainty of 0.2 Hz. 

In the standing wave field, a spatial mismatch of the 
light shift due to (i) the electric-dipole (El) interac¬ 
tion and (ii) the electric-quadrupole (E2) and magnetic- 
dipole (Ml) interactions introduces an atomic-motion- 
dependent light shift 0 , which we refer to as the multi¬ 
polar effect. In addition, the light shift due to the hyper¬ 
polarizability, coupled with the atomic motion in a lattice 
potential, introduces non-linear light shift 0. Conse- 
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quently, the total light shift shows intricate non-linear 
response to the light intensity Ii,. 

We estimate the leading hyperpolarizability shift A/3 /l 
to be 2 X 10“^^ 0 at our maximum lattice laser intensity 
of /l = 89 kW/cm^, which is smaller than our measure¬ 
ment uncertainties of ^ 10“^^. Therefore we neglect the 
other hyperpolarizability effects than the leading term 
and approximate the clock shift as 0, 


hAPc{lL,5Ph,n) 


dAa^^ 

dv-L 


5vi,I-L{r) - A/3JL(r)^ 


+ 


~dAQ^^ 

dv-L 


- Aa'i'" 



-E’r / X 

,E1 


a, 


(1) 


Here, n is the vibrational state for the axial motion, 
and are the differences of El and com¬ 

bined E2-1-M1 polarizabilities of the two clock states, 
and Sui, = I'l, — I'f} is the lattice laser detuning from 
the “El-magic frequency” that equalizes El po¬ 
larizabilities for the ground (g) and excited (e) states, 
i.e., 0 ^ 0 ^ = ^g(e)' note that the Gaussian intensity 

profile /l(^) = of the lattice laser confines 

atoms radially, where wq is the beam radius and r is 
the radial displacement. Eor the typical lattice poten¬ 
tial depth 1 /l ~ 43Er, the radial and the axial vibra¬ 
tional frequencies are ^ 100Hz 

and = yj2Ui,/{mXl) r\j 100 kHz. Since the typical ki¬ 
netic energy ^ h x 100 kHz of the lattice trapped atoms 
is about 10^ times higher than the radial vibrational en¬ 
ergy separation hz^r, we treat the radial motion classi¬ 
cally. In addition, when calculating an energy eigenvalue 
(4 l(^),'^) for the n-th axial vibrational state, the 
axial and radial motion may be decoupled as the adia- 
baticity z^a ^ condition is satisfied. 

To evaluate the light shift given in Eq. (HD , atomic mo¬ 
tion in the lattice plays a crucial role, as the axial motion 
determines the averaged motional state (n) and the radial 
motion determines the effective lattice intensity (/l(^)) 
via the averaged atomic distribution ) in the Gaussian 
intensity profile. The axial vibrational population in the 
n = 0 state is measured to be 0.7(1) by the difference 
of the total area of the red and blue motional sideband 
spectra 0 . As atoms in n > 3 are removed in the 
state preparation, the average occupation is estimated to 
be (n) = 0.4(2) by assuming the atoms are Boltzmann- 
distributed among n = 0,1,2 vibrational states. The 
radial atomic distribution can be inferred from the in- 
homogeneously broadened sideband lineshapes [ 2 ^, as 
the sideband frequency becomes smaller as increases 
due to the Gaussian intensity distribution of 

the lattice laser. The blue-sideband frequency u\^{n,r) = 
[£’J'^(/L(r),n -h 1) - £’g'^(/L(r),n)]/h is given by, 

Ph{n,r) ^ - EB.{n+l)/h. (2) 

Using the sideband spectrum, we determine = 

0 . 8 ( 1 ). 


TABLE I. Corrections and uncertainties for z^Hg/z^Sr measure¬ 
ment. All numbers are shown in a fractional unit of 10“^^. 


Effect 

Correction 

Uncertainty 

Lattice light 

3.6 

6.1 

BBR 

16.1 

3.3 

Atom density 

1.6 

1.6 

2nd-order Zeeman 

6.1 

0.9 

Clock light 

0 

< 0.1 

AOM chirp 

0 

< 1 

Servo error 

-0.4 

0.3 

Hg clock systematic total 

26.9 

7.2 

Sr clock systematic 

-17.0 

0.7 

Gravitational redshift 

-0.5 

0.1 

Ist-order Doppler 

0 

< 2 

Statistic 

- 

3.7 

Total 

9.4 

8.4 


We experimentally measure a^^/h = 

5.5(8) kHz/(kW/cm^) from the axial motional side¬ 
band frequency, where the uncertainty is given by the 
measurement uncertainty of the lattice laser intensity. 
We use calculated value Aa'^^/h = 8.25 mHz/(kW/cm^) 
and A/3/h = —2.5/iHz/(kW/cm^)^ [l^. The en¬ 
tire data are then fitted to Eq. ^ (see Eig. [2]), 
where we employ a multiple regression analysis tak¬ 
ing Jl and z/l as explanatory variables and Auc 
as a response variable. This determines the El- 
magic frequency of z/^^ = 826 855 533(9) MHz and 

~ 1.5(4) X 10“^^/(kW/cm^). 


Eor the measurement of Hg clock frequency, we take 
the lattice frequency to be z/l = — 2 MHz = 

826 855 531MHz so that the light shift Az^c(4l,z^l) 


becomes insensitive, i.e.. 


dll. 


—lop 


0, to 


the variation of the lattice intensity around Gp = 
59(10) kW/cm^. The relevant light shift (black curve) 
and its uncertainty (gray shaded area) are shown in 
Eig. El^b). The lattice light shift is estimated to be 
Auc = —0.04(7) Hz corresponding to the fractional fre¬ 
quency shift of —4(6) X 10“^^. In order to further inves¬ 
tigate the lattice light shift, experimental determinations 
of A/3 and Aa*^™ are crucial. 


As listed in Table [H effects other than the lattice light 
shifts give relatively minor contributions to the system¬ 
atic uncertainty of the Hg clock. The BBR shift is esti¬ 
mated to be —1.6(3) X 10“^^ for the ambient temperature 
of 297(3) K, where the uncertainty mainly comes from 
the theoretical polarizabilities with 10% uncertainty 0 . 
The atom density shift is measured by varying the num¬ 
ber of atoms by 50%. A linear fit to the data points 
infer the collisional shift of —2(2) x 10“^^ for N ^ 1000 
atoms in the lattice, corresponding to the atom density 
X 10^^ cm ^ or an average single-lattice-site occu¬ 
pation of 1.2 atoms. The second-order Zeeman shift is 
investigated by varying Bbias in successive measurement 




















cycles, where the magnetic field is measured through the 
first-order Zeeman shift and the g-factor difference in the 
clock transition [2l[. Fitting to Az/c = ~/5z|Bbias|^ yields 
Pz = 1.6(2) Hz/(mT)^, giving a shift of —6(1) x 10“^^ 
at our operating condition of |Bbias| = 0.21 mT. The 
light shift induced by the clock laser is estimated to be 
< 1 mHz [ill . The servo error is estimated to be smaller 
than the statistical uncertainty of the measurement for 
an averaging time r > 100 s by analyzing the error signal 
of the frequency stabilization to the clock transition. 

Sources of uncertainties for the frequency ratio mea¬ 
surement are also listed in Table HI The Sr clock achieves 
the systematic uncertainty of 7 x 10“^^ [5|. The gravi¬ 
tational redshift is estimated from the height difference 
of Ah = 5(1) cm between the two clocks. The first-order 
Doppler shift between the Hg and the Sr clocks is es¬ 
timated to be less than 2 x 10“^^, which is due to the 
temperature drift in the uncompensated optical path of 
about 5 m. 

Figure [3l^a) displays the Allan standard deviation for 
the ratio measurement z^Hg/^Sr, which shows a cry(r) = 
3 X IO^^^/a/t/s trend for an averaging time of r > 
10 s, mainly responsible to the Dick-effect-limited insta¬ 
bility [i^ of 2 X 10“^^/Y^r/s. The quantum 

projection noise (QPN) 0 only contributes to the in¬ 
stability of « 4 X 10“^^/^r/s. Application of a 

clock laser with smaller instability or synchronous oper¬ 
ation of two clocks to reject frequency noise of the clock 
laser [25| will allow approaching the QPN-limited insta¬ 
bility. 

Figure O^b) summarizes the frequency ratio measure¬ 
ment z^Hg/^Sr- For the evaluation of the statistical uncer¬ 
tainty of z^Hg/z^Sr, we construct a histogram of all the data 
with the bin size of 5 x 10“^^ for a fractional frequency. 
The histogram is fitted to the normal distribution func¬ 
tion, where the reduced chi-squared yields \^Xred = 4.0. 
We thus conservatively evaluate the statistical uncer¬ 
tainty to be 3.7 X 10“^^ by inflating the standard un¬ 
certainty of the mean by \/x^. A weighted mean of the 
result is z/Hg/^Sr = 2.629 314 209 898 909 60(22), where 
the fractional uncertainty of 8.4 x 10“^^ is essentially 
given by the systematic uncertainty of the Hg clock (see 
Table H]). The ratio both of which are 

measured by referencing Cs primary standards at LNE- 
SYRTE [H, EH, is shown by a triangle in Eig. [D^c). This 
result is consistent with our measurements to about la 
uncertainty. Applying the recommended value ^ of 
^^Sr as a secondary representation of the second [8|, the 
absolute frequency of ^^^Hg is given by ^ = 

1 128 575 290 808 155.4(1.1) Hz in the unitTf SI second, 
where an uncertainty of 1 x 10“^^ is given by that of 

,,CIPM 

^Sr 

In summary, we have investigated an optical lattice 
clock based on ^^^Hg and achieved a total systematic 
uncertainty of 7.2 x 10“^^. We determine the magic 
frequency by including the multipolar effect that seri- 



FIG. 3. (a) The Allan standard deviation for the ratio 

measurement z^Hg/z^Sr (circles). Solid line shows a cr^(T) = 
3 X 10“^^/'\/t/s trend. Dashed and dotted lines show the 
Dick-effect-limited and QPN-limited instabilities, (b) The ra¬ 
tio measurements carried out with a 3-month-interval, where 
the solid and dashed lines show the weighted mean and the 
total uncertainty of our measurements. (c) Summary of 
the ratio measurements. Our result (red circle) is displayed 
with that reported by LNE-SYRTE group (blue triangle) 
^SYRTE/^SYRTE indicate Icr uncertainties. 


ously affects the uncertainty budget at low 10“^^. Sim¬ 
ilar strategy will be applied for Sr and Yb-based clocks, 
where the multipolar and hyperpolarizability effects be¬ 
come relevant for the systematic uncertainty at the 10“^^ 
level 0 . We have determined the frequency ratio be¬ 
tween optical lattice clocks z^Hg/^Sr with an uncertainty 
of 8.4 X 10 Accurate determinations of such ra¬ 
tios allow the investigation of the constancy of the fine- 
structure constant a. Taking the Sr clock as an a- 
insensitive anchor, the fractional change Au of the Hg 
clock reveals Aa/a = 0.8Az//z/Hg |l3|. w hich will be 
competitive to the previous constraints [26| if Az^/z^Hg ~ 
10“^^ is investigated over a year. 
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